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We study the process of four-wave mixing (4WM) in ultracold degenerate atomic gases. In particular 
we address the problem of 4WM in boson-fermion mixtures. We develop an approximate description 
of such processes using asymptotic analysis of high order perturbation theory taking into account 
quantum statistics. We perform also numerical simulations of 4WM in boson-fermion mixtures and 
obtain an analytic and numerical estimate of the efficiency of the process. 
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In the recent years atom optics has become a flourish- 
ing subject. With the successful experiments on Bose- 
Einstein condensation JljH new possibilities to study lin- 
ear jjj and nonlinear atom-optics j^j with macroscopic 
wave packets have emerged. So far, the most spectac- 
ular experiment in nonlinear atom optics concerns the 
observation of the so called four- wave mixing of bosonic 
matter waves || . In this nonlinear process, three macro- 
scopic matter wave packets interact and produce a fourth 
one. Also recently, Jin et al. (7) have managed to trap 
and cool a sample of spin polarized 40 K below the Fermi 
temperature in a harmonic trap. This has triggered an 
outburst of activities in trying to understand the proper- 
ties of ultracold fermionic systems Q , and has stimulated 
a great interest in the studies of fermion-boson mixtures 
||. In particular, the question "is non-linear atom optics 
with fermions possible?" has been posed |Tc| ]. 

In this Letter we investigate the process of four-wave 
mixing (4WM) in fermion-boson mixtures and demon- 
strate that under an appropriate choice of parameters 
it is possible to create a macroscopic fourth wave of 
fermions. The recent experiment of Deng et al. with 
Bose-Einstein condensates have demonstrated that 4WM 
in bosonic gases is a truly macroscopic process with an 
efficiency of the order of 6%. In this experiment Bragg 
pulses Q were used to create three condensate clouds, 
which then interacted through collisions. Four wave mix- 
ing of matter waves can also be described - using the 
language of optics - as Bragg scattering from a grating. 
In this picture, two counter propagating condensates cre- 
ate the grating from which the third condensate scatters, 
generating a fourth one. 

An intuitive understanding of Bragg scattering for 
bosons may be obtained considering a homogeneous con- 
densate in a box of volume V. For such a case, the Hamil- 
tonian has the form: 
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is the interaction strength propor- 
tional to the s-wave scattering length a and eg is the 



kinetic energy. Let us consider as an initial state \i) — 
\Ni, ki; JVq, fe; N3, £3) representing Ni particles of mo- 
mentum fej, i = 1,2,3. Assume for the moment that only 
one particle is scattered leading to the final state \ f) — 
\N{,ki; N^,k 2 ] N3,k 3 -,l,ki) with fc 4 ^ h for i = 1,2,3. 
Among all the processes which conserve momentum and 
energy the one corresponding to 



fci 
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is particularly favorable and represents in fact the first or- 
der Bragg scattering |l^,|l3| (see Fig. [IJ) . The above pro- 
cess corresponds to the term at at or a,- in the Hamil- 

1 fc 4 k 2 fe l fe 3 

tonian. This term accounts for the creation of a particle 
in the "grating state" &2 which is already macroscopically 
occupied. This leads to the appearance of a bosonic en- 
hancement factor \/N2 + 1 in the transition amplitude 
associated with the process \i) — > |/). From this point 
of view, due to bosonic enhancement Bragg scattering is 
the most probable process. 

Let us extend the previous analysis to the case 
where many particles are scattered. Clearly, the 
most probable final state is of the form |/) = 
\N[,ki;N2,k 2 ;N^,k3;N4,k4}, where £4 is given by 
Eq.(||). This final state must respect the conservation 
of momentum 



(3) 



and the conservation of the particle number 
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Eqs.(H),(||) and (|j) imply conservation of energy. In the 
reference frame in which k\ and k 2 are collinear the scat- 
tering process is planar. Thus we have five equations 
with six unknown parameters N[, Ni,k±. The actual 
value of 7V4 can be determined by maximizing the tran- 
sition probability Pj/. Using perturbation theory with 
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respect to the off-diagonal elements of the Hamiltonian 
in Eq.[l], we calculate the transition amplitude Tif. For 
Ni = N 2 = N 3 this quantity exhibits a divergence for 
7] = Ni/Ntot = 1/6 w 16.7%, where N tot is the total 
number of particles. This divergence appears due to the 
simplicity of the model and can be removed by using wave 
packets instead of plane waves. In any case, Pif will be 
strongly peaked at N4 ~ N to t/Q- We expect that the 
actual efficiency of the process - which in principle has 
to be calculated non perturbatively - will be of the same 
order of magnitude. Furthermore, our analysis predicts 
a saturation of the efficiency with increasing number of 
atoms, which can be understood as a consequence of the 
Bose statistics. These results are in agreement with the 
experimental results of Ref . || . 

We turn now to the case of fermions. The previous 
analysis clearly stresses the role of bosonic enhancement; 
the macroscopic occupation of the states that form the 
grating select Bragg scattering as the most favorable pro- 
cess. From this point of view, the use of a fermionic grat- 
ing would lead to a rather poor Bragg scattering. For 
this reason we consider Bragg scattering of an incoming 
cloud of fermions on a bosonic grating. In our model, 
fermions and bosons interact via two-body s-wave scat- 
tering. The incoming fermions must obey momentum 
and energy conservation as in the case of boson-boson 
scattering. For fermions there will be inevitably a mo- 
mentum spread due to Fermi statistics. In order to fulfill 
the Bragg condition, illustrated in Fig. [I], the momentum 
spread in the fermionic sample Ak, must be restricted to 
Afc <C I fei — /eg I - This condition can be achieved by trap- 
ping the fermions initially in a sufficiently shallow trap 
which results in a well localized momentum distribution. 

Solving the complete many-body scattering problem 
for the fermions is a formidable task. In our case we can 
simplify the situation assuming first that the fermions 
do not interact. This is a good approximation as long 
as we consider a polarized Fermi gas. The interaction of 
fermions of only one specie is entirely due to p-wave scat- 
tering which is negligible at very low temperatures 
Second, we model the bosonic grating by a scalar poten- 
tial proportional to the local condensate density. In doing 
this we sacrifice in some sense the effects of the statistics. 
In fact, the Fermi statistics appears in our model only 
through the preparation of the initial state, which must 
obey Pauli principle. We also neglect any back-action 
from the fermions on the bosons. This last approxi- 
mation is valid in a situation where the boson density 
is much larger than the fermion density. Finally, since 
the dynamics is basically two dimensional (see Fig. 1) 
we restrict our numerical simulations to 2D. With these 
assumptions we solve the Schrodinger equation for the 
fermions with an effective external potential of the form: 



with ubf 
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Here a denotes the fcrmion-boson 



Vf(x,y) = u B Fn b {x,y) 



(5) 



s-wave scattering length, m the (same) mass for both 
fermions and bosons, L the thickness of the cloud in z- 
direction, whereas ni,(x, y) the equilibrium density of the 
condensate. For the bosons, the trap potential is a com- 
bination of a harmonic potential in the ^direction and a 
periodic structure in x-direction created for instance by 
a standing- wave laser, 
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Vb(x,y) = -m£L c y 



+ Ul cos 2 (|fci — lf2\x). 



(6) 



We have solved numerically the Gross-Pitaevskii equa- 
tion with the potential in Eq. |^, in order to determine 
the equilibrium bosonic density 7^(2;, y), which, in turn, 
results in a periodic potential Vf{x,y) for the fermions. 
The back action of fermions on bosons can be neglected 
provided UBFnf(x,y)/fi <C 1 where (i is the chemical 
potential of the trapped condensate and nf(x,y) the 
fermion density. In our simulation this ratio was kept 
smaller than 0.01. 

Our procedure is now as follows: initially the fermions 
are trapped in a 2D potential of the form V(x, y) — 
mfl 2 (x 2 + y 2 ) centered at (x ,y ). The number of 
fermions is such that the Fermi level still fulfills the con- 
dition fcjr <C \ki — k%\. The trap is then removed and a 
momentum kick hk^ is given to the fermions by illuminat- 
ing the sample with a detuned laser propagating in the 
direction of £3 ||[l5| so the fermion cloud scatters from 
the grating. We finally monitor the density of fermions 
to obtain the efficiency of the process. 

The wave functions of the non interacting fermions ful- 
fill all the same Schrodinger equation 



dt 
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V 2 + V f (x,y) 
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but with different initial conditions for each fermion. The 
initial states are the eigenstates of the displaced harmonic 
potential used for trapping the fermions, 



^n x ,n v (X, y) — A nxjTli 
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-^ 2 ^H nx (x)H n (y), 
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where x = (x — xq)J \Jh/m£l and y — (y — yo)/ y/h/mQ,, 
H nx denotes the Hermite polynomials and A nxtrly is the 
normalization constant. 

The scattering of the fermions can now be numerically 
simulated one by one. We monitor the density in the re- 
gion where the Bragg scattered fermions are supposed to 
appear. In our model, for the fermion cloud we use 40 K 
atoms and a trap frequency of fl/2ir — 10 Hz while for 
the bosonic trap we use a frequency of f2 c = 40O. This 
produces a narrow grating compared to the size of the 
fermionic cloud. Furthermore, we have taken the scat- 
tering length value as a = 6.0 nm and L = 30 /jm with 
N c = 2 • 10 5 bosons in the condensate. The grating has 
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an associate wave number of fci — 1.3 /an -1 (fci = — k 2 ) 
and fei?/fci < 0.7. The initial fermion cloud is positioned 
at xq = 2/0 = ~40 /im and the momentum kick is settled 
to fc 3 = (1.3, 1.3)/im _1 . 

In Fig .|] we show the efficiency r\ = N4/N3 of the Bragg 
process as a function of the total number of fermions 
N3. With an increasing number of particles (N3) the 
efficiency decreases due to the increase of the momen- 
tum spread of the fermions. Fig. || displays a snapshot 
of the cloud after the scattering. We observe that ap- 
proximately 3% of the cloud is scattered in the Bragg 
direction. In this figure one can also see that a part of 
the fermionic cloud is reflected due to the chosen rela- 
tion between the incoming fermion kinetic energy and 
the contrast together with the maximum of the grating 
potential. Note also the appearance of the reflected wave 
packet in the direction — ^3 which in fact corresponds to 
a Bragg reflection. Reflections are drastically reduced if 
the scattering length is negative since in this case the 
potential Vf(x,y) becomes attractive. The Bragg scat- 
tering relies on the periodicity and contrast of the grating 
and it is present for both positive and negative scattering 
lengths. Also, for a fixed interaction time the efficiency 
of the Bragg process decreases with the contrast of the 
grating. 

In order to understand the above results we have de- 
veloped an approach similar to the one used for bosons, 
based on pcrturbativc analysis using plane waves. We 
want to study the efficiency of the process starting with 
the fermionic state \i) = |1, ka + Xii •••> 1> &3 + Xn 3 ) an d 
considering a generic final state in which N4 holes are 
created in the Fermi sea of states with momentum close 
to £3. There are j = (^) different orthogonal final states 
In order to simplify the problem we assume that the 
transition amplitude does not depend strongly on \fj), 
which is a valid assumption as long as Afc <C | fci — fca j . 
We can then evaluate the probability P(N4) of scatter- 
ing N4 particles by taking an average value of the tran- 
sition amplitude so that P(iV 4 ) « (^)|% | 2 - The mean 
efficiency of the process defined as r\ = (N^/N^ is cal- 
culated then using the previous probability distribution 
P(A^). We stress that the aim of this calculation is not to 
obtain an exact expression for the efficiency, but rather 
to estimate whether the scattering process can change 
macroscopically a many particle fermionic state. 

We begin by calculating the transition amplitude 
between the initial state |i) and the final state which 
corresponds to depletion of the first A4 states of the 
Fermi sea around A^. The macroscopically populated 
states of the bosonic grating provide a bosonic enhance- 
ment factor y/ (N 2 + p) [p = 1, 2, N4). Using a time 



dependent perturbative method 16 one obtains 
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where the bosonic enhancement amounts to 
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The energies of the intermediate states (corresponding to 
the scattering of p-fcrmions) are given by 



E p = Eb 



with Ubf 



N 3 



^\ e (£s+xi) 2 +£(*;+xi) 2 
j=p+i j=i 



(ii) 



^y 1 , Ei = E p=0 and k 2 = k 2 e x = -k\e x . 
We denote by Eb the energy of the bosons which is con- 
stant during the process. This Tif 1 is the transition am- 
plitude corresponding to one of the A/4! possible paths 
going from \i) to Hence there are A/4! different Tif 1 's 
and we estimate their contribution by taking an average 
value. To this aim we consider the x ' s as independent 
random variables distributed between —kp and kp with 
a flat probability distribution cx 6{kp — x)- Neglecting 
correlations between different x ' s is a valid approxima- 
tion if 7V4 <C N3. Consequently we consider this regime 
in the following, and restrict ourselves to the values of 
7V4 < O.IN3. The resulting probability of having A'4 
Bragg scattered fermions is 



P(N 4 ) cx ( -1 
or 



N3—N4 



1 



(A^3 — A'4)! 



(12) 



where we have assumed A^i = N% . This Poissonian prob- 
ability distribution depends on the number of fermions 
A^3 and on the parameter 



(13) 
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which is also A^-dependent through the dependence on 
the Fermi momentum kp. The parameter a contains also 
the boson density rife. In Fig. || we display the efficiency 
r), (defined as the mean value of A4 over A3) as a function 

1 /3 

of N3 and the A^-independent parameter aN 3 ' . We 
observe that for small values of a the efficiency decreases 
with increasing A3. On the contrary, with larger a the 
efficiency can reach larger values, only limited by the 
assumptions of the model, i.e rj < 0.1. 

In summary, we have discussed the effects of Bose and 
Fermi statistics on four- wave mixing processes. For a 
pure bosonic process the Bose statistics sets a funda- 
mental limit in the efficiency. In the case of an incoming 
fermionic cloud, the result depends strongly on the val- 
ues of the various physical parameters involved in the 
problem. On one hand, we have performed a numer- 
ical analysis which shows that the momentum spread 
A A; <^ I fci — k 2 \ is a crucial parameter in order to ob- 
tain a macroscopic fourth wave. On the other hand, the 
analytical treatment which assumes that Bragg condi- 
tion is always fulfilled, exhibits an interplay between the 



3 



statistical and collisional effects leading to a decreasing 
efficiency for small N$ and increasing efficiency for large 
N3. Typically, for small a, the value of {N4) is negligible 
meaning a small efficiency On the contrary, for suffi- 
ciently big values of a, the efficiency is only limited by 
the assumptions of the model, namely the lack of correla- 
tions between the scattered particles. The efficiency can 
then reach values of the order of a few percent which char- 
acterizes the creation of a macroscopic fermionic fourth 
wave. 

The authors wish to thank W. Ketterle for helpful 
discussions. We acknowledge support from Deutsche 
Forschungsgemeinschaft (SFB 407) and the TMR net- 
work ERBXTCT96-0002. 
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FIG. 2. Numerical estimate of the efficiency 77 of the 4WM 
versus the number of incoming fermions N3. r\ decreases for 
larger values of N3 due to the spread of the momentum for 
the incoming fermions. 



[1] 
[2] 
[3] 
[4] 
[5] 

[6] 
[7] 
[8] 



[9] 
[10] 

[11] 
[12] 

[MI 



M.H. Anderson et ai, Science 269, 198 (1995) 

K.B. Davis et ai, Phys. Rev. Lett., 75, 3969 (1995) 

C.C. Bradley et al, Phys. Rev. Lett. 75, 1687 (1995) 

K. Bongs et al, Phys. Rev. Lett. 83, 3577 (1999) 

G. Lenz, P. Meystre, and E.W. Wright, Phys. Rev. Lett. 

71, 3271 (1993) 

L. Deng et al, Nature 398, 218 (1999) 

B. DeMarco and D.S. Jin, S cience 285, 1703 (1 999) 



M.A. Baranov, D. S. Petrov. 



cond-mat/9901108 



L. Vichi 



and S. Stringari, Phys. Rev. A, 60, 4734 (1999), H.T.C. 
Stoof and M. Houbiers, Proceedings of the International 
School of Physics "Enrico Fermi" Course CXL, eds. M. 
Inguscio, S. Stringari, and C.E. Wieman (IOP press Am- 
sterdam 1999). 

N. Nygaard, K. Moelmer, Phys. Rev. A, 59, 2974 (1999) 
W. Phillips, private communication. 
M. Kozuma et al, Phys. Rev. Lett. 82, 871 (1999) 
M. Trippenbach, Y.B. Band, and P.S. Julienne, Optics 
Express 3, 530 (1998) 
M. Trippenba ch, Y.B. 
mat/0002119 



Band, and P.S. Julienne, cond 



[14] 
[15] 
[16] 



B. DeMarco et ai, Phys. Rev. Lett. 82, 4208 (1999) 
S. Inouye et ai, Science 285, 571 (1999) 
LD. Landau and E.M. Lifshitz, Quantum Mechanics 
(Pergamon Press, London-Paris, 1958) 




kj -k 2 

FIG. 1. In the reference frame in which fci and are 
collinear the scattering process can always be described in 
a plane. The criterion £4 = fci — £2 + £3 together with the 
conservation of energy \k^\ = \k?,\ gives the Bragg condition 
kix = ki with fci = — &2 chosen in the direction of the x-axis. 
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FIG. 3. A Snapshot of the density of the fermionic cloud 
after the scattering. In this case 3% of the atoms are found 
in the Bragg direction. The grating is represented by the two 
counter propagating waves fci = —k^ 
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FIG. 4. Analytical estimate of the Bragg process efficiency. 
The model is limited by the assumption 77 <C 1. For small 
values of a the efficiency decreases with increasing N3. 
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